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Feedback  constant  — front  panel  setting 

Internal  feedback  constant 

ANC  filter  length 
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Signal  (sinusoid)  RMS  power 
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I.  INTRODUCTION 

This  report  describes  tests  performed  on  the  Adaptive  Noise  Canceller  (ANC)  as 
part  of  a joint  NUSC/NOSC  program  investigating  adaptive  processors.  The  purpose  of  the 
test  was  to  evaluate  the  effects  of  various  inputs  and  internal  parameters  on  the  operating 
characteristics  of  the  ANC  (fig  1 ). 

The  primary  dependent  parameters  of  interest  were  the  adaptive  time  constant  r, 
the  notch  width  W,  and  the  notch  depth  D.  These  were  measured  as  a function  of  the  feed- 
back constant  p,  the  filter  length  N,  the  sample  rate  f^,  the  reference  power 
the  delay  A.  In  one  experiment,  the  reference  signal-to-noise  ratio  was  varied  as  well.  The 
testing  also  included  measurement  of  the  ANC  transfer  function  and  the  input-output 
cross-correlation. 

According  to  theory  (ref  1),  the  shortest  time  constant  of  the  ANC  for  cancelling 
sine  waves  in  white  noise  is  given  by 

T=[2^(-y  a2  + 02)f^]-l;  d) 

where, 

T = Decay  time  constant  (sec)  for  the  output  in  volts  RMS  (this  is  twice 
the  decay  time  for  the  output  power), 

Sample  rate. 

Feedback  constant  in  (Volts)~2^ 

RMS  power  in  (Volts)^  for  the  narrowband  components  in  the  reference 
input, 

RMS  power  of  the  broadband  component  of  the  reference  input. 

Special  caution  must  be  taken  to  be  sure  that  the  bandwidth  of  the  “white”  noise  is  exactly 

Nyquist;  i.e.,  fg/2  (see  appendix  A).  • , 

It  should  also  be  noted  that  the  above  theory  assumes  convergence  of  the  weights  | 

to  a linear  filter.  For  a sinusoidal  input  this  is  not  always  the  case  (ref  2);  the  solution  for  [ 

the  weights  may  not  be  stationary  and  one  might  even  question  the  definition  of  “con-  i 

vergence.”  However,  such  problems  are  beyond  the  scope  of  this  study. 


1.  Treichler  JR,  The  Spectral  Line  Enhance,  Ph.D  Dissertation,  Stanford  University,  Stanford,  California, 
May  1977 

2.  (NOSC)  NUC  TN  1617,  Adaptive  Noise  Cancelling  of  Sinusoidal  Interference  (U),  by  J Glover,  Unclas- 
sified, December  1975 
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For  a narrowband  reference  signal  {a}  = 0),  the  notch  width  W,  is  given  (ref  2) 


W = ja  N /tt, 


(2) 


where  N is  the  filter  length.  Equations  ( 1 ) and  (2)  imply 


W = — . 

TTT 


(3) 


The  input  to  the  primary  and  hence  the  delay,  should  have  no  effect  on  t for  narrowband 
signals  (ref  2).  However,  in  the  broadband  case,  the  finite  filter  length  must  be  taken  into 
account  (ref  3).  In  fact,  for  Ap  > N there  will  be  almost  no  cancellation.  The  optimal 
value  of  Ap  is  about  N/2  (ref  3). 

Four  tests  were  performed.  The  first  examined  broadband  cancellation  (ff|  = 0). 
The  time  constant  t was  measured  as  a function  of  n (test  1.1),  N (test  1.2),  Ap  (test  1 .3), 
and  (test  1.4)  and  then  compared  with  equation  (1). 

The  second  test  created  narrowband  cancellation  (o^  “ measuring  the  time 
constant  as  a function  of  m (test  2.1),  N (test  2.2),  and  (test  2.3).  The  notch  widths 
were  also  measured  as  a function  of  fi  and  N in  experiments  2.4  (signal  plus  noise  in 
primary)  and  2.5  (noise  alone  in  primary). 

In  the  third  test,  multiple  lines  were  present.  Test  3.1  contained  two  sinusoids  in 
the  reference;  test  3.2  had  two  sinusoids  in  the  primary;  and,  test  3.3  had  four  sinusoids  in 
both  the  primary  and  reference.  Test  3.4  was  similar  to  3.3;  however,  the  sinusoids  were 
input  to  four  distinct  references. 

Test  Series  4 investigated  the  effects  of  adding  an  independent  noise  source  to  the 
reference.  The  first  test,  4.1,  measured  the  time  constant  for  a reference  consisting  of  a 
sinusoid  plus  white  noise.  The  reference  signal-to-noise  ratio  was  varied  from  oo  to  -18  dB. 
Test  4.2  measured  notch  width  and  notch  depth  as  a function  of  reference  SNR  with  both 
rsiial  and  noise  in  the  primary.  Test  4.3  repeated  test  4.2  with  noise  alone  in  tlie  primary. 


U.  TEST  SERIES  1 - BROADBAND  CANCELLATION 


TEST  OBJECTIVE 

The  purpose  of  this  series  of  tests  was  to  examine  the  performance  of  the  ANC  for 
primary  and  reference  inputs  of  identical  broadband  noise.  In  parti  :ular,  the  adaptive  time 
constant  was  measured  as  a function  of  pi,  N,  Ap,  and  and  then  compared  to  that 

predicted  by  theory  (equations  ( 1 ) and  (A-10)). 


3.  Widrow  B,  Glover  J,  McCooi  J,  et  al.  Adaptive  Noise  Cancelling:  Principles  and  Applications,  Proc  IEEE 
vol  63,  December  1975 
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TEST  DESCRIPTION 


The  experimental  setup  for  this  series  of  tests  is  pictured  in  figure  2.  The  input  to 
the  primary  and  reference  consisted  of  the  same  white  Gaussian  noise,  low-pass  filtered  at 
600  Hz.  Prior  to  each  run,  the  ANC  filter  weights  were  initialized  to  zero.  A spectrum 
analyzer  and  X-Y  plotter  were  used  to  plot  the  input  (primary)  and  output  spectra  both 
before  and  after  cancellation. 

During  convergence,  the  output  voltage  was  recorded  with  an  RMS  voltmeter  and 
plotted  as  a function  of  time.  The  time  constant  t was  measured  by  fitting  the  curve 
Ae~*/^  + Aq  to  this  cu'put.  Since  Ag,  the  minimum  mean  square  error,  was  small,  this 
was  approximately  the  time  required  for  the  ANC  output  to  drop  to  (1/e)  of  its  initial 
value  (see  fig  3). 

A summary  of  the  parameter  values  used  in  the  various  tests  in  this  series  appear 
in  figure  4.  For  each  test,  a plot  was  made  of  r versus  the  pertinent  dependent  parameter 
and  compared  to  equation  (1). 


TEST  RESULTS 

Test  1 . 1 measured  r as  a function  of  n for  fixed  N = 64,  Ap  = 32,  and  = 1 .56. 

The  theoretical  dependence  is  given  in  appendix  A,  equation  (A-10).  Figure  5,  a plot  of 
the  results,  shows  excellent  agreement  between  experiment  and  theory. 

In  test  1.2,  T was  measured  as  a function  of  N with  fixed  n and  was  set 

equal  to  N/2.  The  reasons  for  this  are  discussed  in  the  following  paragraph.  As  seen  from 
equation  (A-10),  r should  exhibit  no  dependence  on  N.  This  was  borne  out  by  the  results, 
which  are  plotted  in  figure  6. 

Test  1 .3  examined  the  performance  of  the  ANC  as  a function  of  the  delay  Ap.  The 
other  parameters  were  fixed:  /a  = 2“*  2,  N = 64,  = 1 .58.  With  an  infinite  filter  length, 

the  ANC  would  show  no  dependence  on  Ap;  however,  a finite-length  filter  cannot  model 
arbitrarily  long  delays  (ref  3).  Figure  7 contains  a plot  of  the  magnitude  of  the  ANC  trans- 
fer function  at  500  Hz  versus  the  delay.  Because  of  end  effects,  one  would  expect  optimal 
cancellation  for  Ap  = N/2,  one  half  the  filter  length.  For  Ap  > N,  the  delay  cannot  be 
modelled  by  the  filter  and  no  cancellation  occurs.  On  the  other  hand,  Ap  should  not  affect 
the  time  constant  of  the  ANC  provided  it  is  measurable.  This  is  possible  as  long  as  there 
is  some  observed  cancellation.  These  remarks  are  in  agreement  with  figure  8,  which  plots 
T versus  A. 

In  test  1.4,  T was  measured  as  a function  of  the  reference  voltage  The 

parameters  ju,  N,  and  Ap  were  fixed  at  2"'  64,  and  32  respectively.  The  reference  voltage 

was  incremented  at  steps  of  3 dB.  For  values  below  -1 2 dB,  the  cancellation  was  incom- 
plete. This  was  caused  by  the  inability  of  the  filter  weights  to  provide  enough  gain  to  raise 
the  reference  input  to  the  level  of  the  primary.  For  these  values,  it  was  not  possible  to 
accurately  measure  t.  Tlie  results  are  plotted  in  figure  9 and  once  again  show  excellent 
agreement  with  theory. 


AMPLITUDE 


PRIMARY 


Figure  2.  Experimental  setup  for  Test  Series  1 . 


i 


9 


Test  1.  Broadband  Noise  Cancellation 


1.1  T versus  fi 


N 


= 2~®  2”^  2”^®  2*^^  2*^^  2”^^  2”^^ 


= 64 
= 32 


'^RMS  = 1-56  volts 


1.2  T versus  N 


= 2-12 


= 16,64,256,1024,1964 


= N/2 


'^RMS  = 1-58  volts 


1.3  T versus  A 


N 


= 2-12 


= 64 


Ap  = 0,  16,  32,  63,  64,  65,  128 


'RMS 


= 1.58  volts 


1 .4  T versus  V 


RMS 


U 

N 

A„ 


= 2-12 


= 64 
= 32 


Vrms  ° -95,  .66,  .47,  .33,  .23,  .17,  .12,  .087,  .064,  .049,  .038,  .031 


Figure  4.  Oulline  of  parameters  for  test  series  J.  The  sample  frequency  f was  fixed  at  2048  Hz. 
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Figure  5.  Test  1.1  — plot  of  r versus  fi. 


» ■ ■ MEASURED 

THEORY  ^ 


Figure  7.  Test  1.3  — plot  of  ANC  transfer  function  magnitude  at  500 
Hz  versus  delay  (in  units  of  filter  taps)  for  N = 64. 

— MEASURED  , . ij 


THEORY 
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III.  TEST  SERIES  2 - NARROWBAND  CANCELLATION 


TEST  OBJECTIVE 

The  second  series  of  tests  was  designed  to  examine  the  performance  of  the  ANC 
when  the  reference  signal  consists  of  a single  sinusoid.  The  adaptive  time  constant  t and 
the  notch  width  W were  measured  for  various  values  of  n,  N,  and  in  an  attempt  to 

verify  equations  ( 1 ),  (2),  and  (3). 


TEST  DESCRIPTION 

The  experimental  setup  for  this  series  of  tests  is  pictured  in  figure  10.  The  measure- 
ment of  the  time  constant  t was  identical  to  that  of  Test  Series  1.  The  notch  widths  were 
measured  at  the  3 dB  down  points  of  the  output  spectrum  (see  figure  1 1 ).  A summary  of 
the  parameter  values  used  in  these  tests  appears  in  figure  1 2. 


TEST  RESULTS 

Tests  2. 1,  2.2,  and  2.3  measured  the  dependence  of  the  time  constant  r on  N, 
and  respectively.  The  relationship  predicted  by  the  theory  is  equation  (A-2)  of 

appendix  A.  The  test  results,  along  with  the  theoretical  curves,  appear  in  figures  13  through 
1 5 and  show  excellent  agreement. 

Test  series  2.4  and  2.5  measured  the  dependence  of  the  notch  width,  W,  on  n and  N. 
In  2.4  the  primary  contained  the  same  sinusoidal  signal  as  the  reference,  plus  additive  noise 
at  a signal-to-noise  ratio  (RMS  power  across  the  entire  ANC  input  band)  of  0 dB.  In  2.5, 
the  primary  contained  only  the  signal. 

Test  2.4. 1 examined  notch  width  as  a function  of  /i.  The  experiment  was  repeated 
for  two  values  of  N:  N = 32  and  N = 5 1 2.  Results  are  plotted  in  figures  16  and  1 7 along 
with  the  theoretical  curves  (equation  A- 13).  In  test  2.4.2,  p was  fixed  at  2“®  and  N was 
varied.  The  large  value  of  p was  chosen  in  order  to  obtain  notches  sufficiently  broad  to 
measure.  These  results  appear  in  figure  18.  Finally,  the  data  points  obtained  are  plotted 
versus  rt^eoretical  figure  19  along  with  equation  (3).  Again,  theory  and  experiment  are 
consistently  in  close  agreement. 

Tests  2.5.1  and  2.5.2  were  identical  to  2.4.1  and  2.4.2,  except  that  the  signal  was 
absent  in  the  primary.  The  results  are  plotted  in  figures  20  through  23.  As  may  be 
observed,  the  absence  of  the  signal  in  the  primary  had  no  effect  on  the  notch  width. 
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dBV  per  0.125  Hz 


Figure  1 1.  Example  of  the  measurement  of  notch  width  (3  dB  down  point)  for  test  senes  2. 
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W versus  N and  p;  primary  signal-to-noise  ratio  = -<» 

2.5.1 

W versus  fi  (two  values  of  N) 
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Identical  to  2.4.1  except  • 1.1  volte 

2.5.2 

W versus  N 

i 

) 

1 

Identical  to  2.4.2  except  =1.8  volte 

Figure  12.  Parameters  for  test  series  2.  The  sample  frequency  was  set  at  fj  = 2048,  the  delay  at  A_  = N/2, 
one-half  the  filter  length  and  the  signal  frequency  225  Hz. 
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Figure  19.  Test  2.4  — plot  of  notch  width  versus  t.  The  data  points  were  taken  from  2.4.1  and 
2.4.2.  For  clarity  a portion  of  the  plot  has  been  redrawn  at  a magnified  scale. 


Figure  20.  Test  2.5.1  - plot  of  notch  width  versus  /z  for  N = 32. 
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Figure  21.  Test  2.5.1  - plot  of 


Figure  22.  Test  2.5.2  - plot  of 
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IV.  TEST  SERIES  3 - MULTIPLE  LINES 


TEST  OBJECTIVE 

This  series  of  tests  examined  the  behavior  of  the  noise  canceller  foi  ; i.Tials  contain- 
ing more  than  one  sinusoid.  Time  constants,  notch  widths,  and  notch  depths  were  measured 
as  a function  of  the  frequency  separation  Af,  of  the  signal  components. 


TEST  DESCRIPTION 


Tlie  experimental  setups  for  Test  Series  3 are  pictured  in  figures  24  through  27. 
Figure  28  is  a diagram  of  the  ANC  with  multiple  references,  which  was  used  in  experiment 
3.4  in  place  of  the  configuration  depicted  in  figure  1 . 


V 

Slightly  different  techniques  were  used  to  measure  the  time  constants  in  this  series 
of  tests.  The  noise  in  the  primary  (which  was  not  cancelled)  overshadowed  the  output 
power  in  test  3.1 . In  order  to  make  t visible,  the  level  of  the  frequency  cell  containing  the 
signal  ( 1 20  Hz)  was  monitored.  The  level  was  in  decibels  and  t was  computed  from  the 


formula  (see  appendix  B)  r ~ 3-3  and  3.4,  monitoring  was  not  neces- 

sary; however,  the  noise  did  tend  to  obscure  the  final  level.  As  a result,  three  points  were 
needed  to  fit  the  curve  + Ag  and  thus  determine  t. 


Notch  widths,  W,  were  measured  as  in  test  series  2.  The  notch  depths,  D,  o(  the 
iranst^  functions  were  (Ipte^ming^  by.com|iari^  the  inpm  ai^  ouj^ut  spectra  of  the 
primary,  (To  remove  the  effects  of  the^xternal  bOtHlz  t%ir7*the  input  spectrum  was 
measured  at  the  output  by  setting  the  reference  voltage  to  zero.)  More  precisely,  the  notch 
depths  are  given  by 


D = primary  output  power 
in  dBV  at  signal  frequency 


primary  input  power 
in  dB  at  signal  frequency. 


The  parameter  values  used  in  Test  Series  3 appear  in  figure  29. 


TEST  RESULTS 

Test  3.1  measured  the  time  constant  t,  notch  width  W,  and  notch  depth  D as  a 
function  of  frequency  separation  Af  of  two  sinusoids  in  the  reference.  One  frequency  was 
held  constant  at  120  Hz  and  the  other  varied  from  170  Hz  to  120.5  Hz.  The  120  Hz  tone 
was  also  present  in  the  primary. 

In  all  cases  two  notches  were  visible,  at  the  two  reference  frequencies.  For  separa- 
tion of  less  than  10  Hz,  the  notch  widths  were  not  measurable.  These  results  are  tabulated 
in  table  1.  The  absolute  notch  depths,  dj,  of  the  primary  output  spectrum  are  included  as 
well  as  the  notch  depths,  Dj,  of  the  transfer  function.  The  depths  d]  and  apPear  to  be 
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Figure  24.  Experimental  setup  for  test  3.1. 
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Figure  25.  Exp'.-rimeiital  setup  for  test  3.2. 


•e  26.  Experimental  setup  for  test  3.3. 
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Figure  27.  Experimental  setup  for  test  3.4. 


of  multiple  reference  ANC. 


TEST  3 


3.1  r,  W,  D versus  Af 


^ = 2-® 

N = 64 

Primary;  = .72  = .75 

RefererKe:  = .78  V\9  = ,78 

f,  = 120  Hz  f2  = 120.5,  121,  122,  123,  125,  130,  140,  150, 170  Hz 


3.2  W versus  Af 


N = 64 

Primary:  = .79  * .39  = .79 

f,  = 120  Hz  = 119.5,  119.75,  120.25,  120.5,  121,  122,  123,  125, 
130,  140,  150,  170  Hz 

Reference:  « .39 


3.3  r,  W,  D for  four  signals  in  reference 


N = 64 

Primary:  = .49  V*;9  = -37  = .36  Vf  = .36  Vj9  = .37 

f,  = 201  Hz  fj  = 201.1  Hz  fg  = 201.2  Hz  f^  = 201.3  Hz 
Reference;  = .36  V^'9  = .75  V^'9  = .51  V|'9  = .41 

3.4  Identical  to  3.3  except  the  four  signals  were  input  to  four  different  references  and  N = 16. 


Figure  29.  Parameters  for  test  series  3.  The  sample  rate  was  set  at  = 2048  and  the  delay  at  A^  = N/2, 
one  half  the  filter  length.  ^ 
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approximately  equal  at  the  two  reference  frequencies.  However,  the  1 20  Hz  signal  in  the 
primary  was  about  36  dB  above  the  noise  level  (at  a resolution  of  .125  Hz).  Thus,  con- 
siderably more  cancellation  occurred  for  the  sinusoid  which  appeared  in  both  the  primary 
and  reference  than  for  the  cancellation  which  occurred  in  the  reference  alone;  i.e.,  I D j I 
was  generally  36  dB  larger  than  I D2 1 . Also  notice  the  decrease  in  cancellation  (at  botii 
f|  and  f2)  with  decreasing  Af  (especially  for  Af  < 10  Hz);  therefore,  the  above  observations 
may  not  be  explained  simply  as  an  infinitely  deep  notch  which  was  clipped  by  the  electronic 
equipment.  I 


Table  I . Notch  Widths  and  Notch  Depths  for  Experiment  3.1 . The  First  Signal  was  at  1 20  Hz 
with  a Power  of  -3  dBV.  The  Noise  Level  was  -38.5  dB/Cell  at  a Cell  Width  of  .125  Hz. 
d Indicates  Absolute  Depth  of  the  Notch;  D Shows  Cancellation  in  dB. 
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Af 

h 

Wl 

(Hz) 

W2 

(Hz) 

dl 

power/cell 

dBV 

d2 

power/cell 

dBV 

Dl 

dBV 

D2 

dBV 

50 

170 

2.8 

— 

-57 

_ 

-54 

_ 

30 

150 

2.0 

2.0 

-59 

-63 

-56 

-25 

20 

140 

1.3 

1.6 

-55 

-56 

-52 

-18 

10 

130 

.8 

1 

-48 

-48 

-45 

-10 

5 

125 

- 

- 

-43 

-41 

-40 

-3 

3 
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- 

- 

-42 

■A2 

-39 

■4 

2 

122 
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-44 

-45 

^1 

-7 

121 

- 

— 1 

-38 

-38 

-35 

0 

.5 

120.5 

NOT  CONVERGED 

According  to  the  theory  in  reference  4,  two  time  constants  should  be  present,  with 
values  given  by 


T 1 

* 2ttaj{^\± 


where 


X^  2 - 2 


and 


Z = 


sin  -f — Af 
‘s 

sin  Y Af 
‘s 


(4) 


4.  NOSC  Working  Paper,  Time  Constants  and  Learning  Curves  of  LMS  Adaptive  Filters  (U),  by  M Shensa, 
Unclassified,  1979 
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For  large  frequency  separations,  the  time  constants  and  t_  are  approximately  equal 
with  values  given  by  substituting  X = N/2.  For  small  Af,  however,  they  are  very  different; 
the  fast  time  constant  has  X = N (twice  as  fast  as  when  Af  is  large).  Note  that  this  is  the 
same  as  (A-2)  with  replaced  by  the  total  signal  power.  The  slow  time  constant  is 

proportional  to  (more  precisely  X_  is  about  7r^N^(A0^/(  1 2 fj )).  j 

1 

Both  time  constants  were  observed  in  the  present  test.  An  example  appears  in 
figure  30.  The  system  response  time  prevented  accurate  measurement  of  the  fast  time 

constant,  however,  considerable  success  was  achieved  for  the  slower  time  constant,  t_.  J 

The  results  appear  in  table  2 and  figure  3 1 (note  that  ^ = p/2).  , 


Table  2.  Time  Constants  Versus  Frequency  Separation  for  Test  3.1 , Theoretical  and  Measured  Values 


Af 

Hz 

u 

sec 

T’. 

sec 

(measured)  sec 

.5 

.08 

401 

530 

1 

.08 

97 

100 

2 

.08 

24 

21 

3 

.088 

6.1 

3 

5 

.083 

2.7 

— 

10 

.09 

1.0 

- 

20 

.11 

.29 

.2 

30 

.13 

.20 

.2 

50 

.13 

.20 

.17 

The  configuration  in  test  3.2  was  essentially  the  same  as  in  test  series  2,  except  for 
the  addition  of  a second  sinusoid  in  the  primary.  Notch  width  was  measured  as  a function 
of  Af,  the  signal  separation  in  the  primary.  In  theory,  the  second  signal  should  have  no 
effect  on  the  notch  width.  This  was  confirmed  by  the  results  (fig  32),  although  it 
should  be  noted  that  the  measurement  accuracy  was  limited  by  the  narrowness  of  the 
notches. 

Test  3.3  examined  the  behavior  of  the  noise  canceller  for  a reference  signal  consist- 
ing of  a sum  of  four  sinusoids  spaced  Af  = . 1 Hz  apart  and  with  relative  signal  powers  of  1 , 

4,  2,  and  1 .33.  Extrapolating  the  theory  from  the  case  of  two  closely  spaced  sinusoids 
(ref  4),  one  would  expect  the  fastest  time  constant  to  be  given  by  t = (p  Nf^)"*  where 
is  the  total  reference  power  and  p =-^,  one  half  the  panel  setting  for  p.  The  actual 
measured  time  constant  for  N = 64  and  p = 2“*2  ^as  1 .8  sec.  The  total  reference  voltage 
was  1.06  volts  RMS.  When  these  values  (o^  = V“/25)  are  substituted  in  the  above  formula, 
one  gets  1 .6  seconds,  in  close  agreement  with  the  measured*  value. 

On  the  other  hand,  the  slow  time  constant  should  be  longer  than  the  fast  time  constant 
/ fs  V 5 

by  a factor  of  1 , which  is  about  10-^  X (1.2)  seconds  or  1.4  days.  For  this 
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reason,  the  slower  time  constant  was  not  observecJ.  However,  its  existence  is  indicated  by 
the  results.  Figure  33  shows  the  supposedly  converged  output  (several  minutes).  The 
cancellation  is  considerably  less  (27  dB  less)  than  in  Test  3.4.  It  was  not  possible  to  meas- 
ure the  notch  width.  The  transfer  function  notch  depth  was  of  the  order  of  -3 1 dB 
- (-20  log  .37  volts)  = -22  dB. 


Test  3.4  was  identical  to  3.3,  except  that  four  separate  reference  inputs  were  used 
for  the  four  sinusoids.  In  this  case,  one  might  guess  that  the  time  constants  wduld  act 
independently.  The  reference  voltages  ranged  from  .36  volts  to  .71  volts.  Since  N was  16 
in  this  test,  the  above  values  give  a range  in  time  constants  (equation  A-2)of  6 seconds 
to  48  seconds.  The  average  power /reference  input  was  .5  volts,  which  corresponds  to  a 
time  constant  of  24  seconds.  The  actual  measured  r was  17  seconds  and  falls  in  the  above 
range  of  values. 

Unlike  test  3.3,  test  3.4  did  not  appear  to  have  a “hidden”  slow  time  constant.  The 
spectrum  of  the  output  (found  in  fig  34)  shows  cancellation  down  to  the  level  of  the 
noise  (as  was  the  case  for  the  converged  output  in  test  3.1 ),  thus  indicating  convergence. 

The  notch  width  could  not  be  measured.  Cancellation  was  of  the  order  of  -58  - (20  log  .32). 


V.  TEST  SERIES  4 - REFERENCE  CONTAMINATED  BY  NOISE 


TEST  OBJECTIVE 

This  series  of  tests  was  designed  to  examine  the  behavior  of  the  noise  canceller  for 
a reference  consisting  of  a single  sinusoid,  plus  noise  independent  of  the  primary.  Time 
constants  t,  notch  widths  W,  and  notch  depths  D,  were  measured  as  a function  of  the  ref- 
erence signal-to-noise  ratio,  SNR. 


TEST  DESCRIPTION 

The  experimental  setup  for  test  series  4 appears  in  figure  35.  The  reference  signal- 
to-noise  ratio,  SNR,  was  measured  as  in  figure  1 , after  the  600  Hz  filter.  The  measurement 
of  the  time  constants  and  notch  widths  followed  the  procedure  of  test  series  2,  while  the 
notch  depths  were  determined  as  in  test  series  3.  The  parameters  used  in  this  series  of  tests 
appear  in  figure  36. 
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Figure  35.  Experimental  setup  for  test  series  4. 
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4.3  Identical  to  4.2  except  for  absence  of  the  signal  in  the  primary. 

Figure  36.  Parameters  for  test  series  4.  The  sample  rate  was  = 2048  and  the  delay,  = N/2,  was  one 
half  the  filter  length. 

TEST  RESULTS 

Test  4.1  measured  the  time  constant  t as  a function  of  reference  signal-to-noise 
ratio,  SNR.  The  results  are  plotted  in  figure  37,  They  are  compared  with  equation  (A-8) 
in  table  3.  Although  the  general  tendency  (longer  time  constants  for  decreasing  SNR)  is 
correct,  the  values  differ  widely  from  equation  (A-8)  for  low  SNR.  This  may  be  due  to  an  ina- 
bility to  measure  the  effects  of  the  noise  on  the  time  constant  (since  the  weights  start  at  zero), 
which  is  the  converged  final  value  for  uncorrelated  noise  alone  in  the  reference.  Or,  it  may  be 
a failure  of  the  theory  to  account  for  the  type  of  behavior  described  in  reference  2. 


Test  4.2  measured  notch  width  and  depth  as  a function  of  signal-to-noise  ratio.  N 
was  fixed  at  64  and  /i  at  2 ~®.  The  primary  contained  both  the  signal  and  noise.  Wiener 
filter  theory  (ref  5)  indicates  that  the  formula  for  notch  width,  equation  (A-13),  should 
remain  unchanged;  whereas,  the  notch  depth  (theoretically /infinite  for  no  noise  in  the 
reference)  should  be  given  by 


D = -20  log  ( 1 


(5) 
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5 . (NOSC)  NUC  Working  Paper,  Adaptive  Noise  Cancelling  and  Line  Enhancing  in  Passive  Sonar  Signal 
Processing:  Theory , Experiment,  and  Hardware  Requirements  (U),  by  P Reeves  and  D Chabries, 
Confidential,  1976 
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Figure  37.  Test  4.1  — plot  of  t versus  SNR,  the  reference  signal-to-noise  ratio. 


Table  3.  t vs  SNR  for  Test  4.1  Compared  to  Values  Computed  by  Equation  (A-8) 


Reference 

T (equation  A-8) 

r (measured) 

SNR  in  dB 

seconds 

seconds 

oo 

1.6 

1.5 

24 

1.6 

1.5 

18 

1.6 

1.4 

12 

1.6 

1.5 

6 

1.6 

1.3 

0 

1.6 

1.4 

-3 

2.7 

-6 

5.8 

2.1 

-9 

11 

2.8 

-12 

19 

2.3 

-18 

48 

12 

10 

The  data  obtained,  plotted  in  tlgure.s  38  and  39,  are  consistent  with  these  remarks. 

Test  4.3  was  identical  to  te.st  4.4  except  that  the  signal  was  absent  in  the  primary.  The 
results  are  plotted  in  figures  40  and  41 . Although  the  notch  width  is  the  same  as  in  Test 
3.2,  the  magnitude  of  the  notch  depth  is  considerably  smaller  (although  still  linear  on  a 
log-log  scale).  This  is  very  similar  to  the  behavior  observed  for  the  two  notches  in  Test  3. 1 . 
The  inability  of  present  theory  to  explain  these  results  points  out  a need  for  additional 
theoretical  analysis. 
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Figure  38.  Test  4.2  - plot  of  notch  width  versus  reference  signai-to-noise  ratio. 


Figure  39.  Test  4.2  - plot  of  notch  depth  versus  reference  signal-to-noise  ratio. 
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Discussed  her^are  some  of  the  interrelations  between  equations  ( 1 ) through  (3)  and 
the  actual  experimental  setup.  The  internal  scaling  of  the  ANC  box  (see  fig  1)  results  in  the 
relationships 


A 2 - 
M =— 


(A-1) 


where  p is  the  feedback  constant  in  (volt)“^,  n is  the  feedback  as  read  on  the  front  panel, 
and  is  the  reference  power  in  (volt)2. 


For  a single  sinusoid,  with  the  weights  initially  set  at  0,  the  output  will  exhibit  a 
single  time  constant  (refs  1, 3,  4)  given  by  the  substhution  of  (A-1)  into  (1) 

(with  = 0) 


T = 


(A-2) 


Note,  that  the  sample  rate  fj,  is  included  to  convert  from  time  in  number  of  sample  points 
to  time  in  seconds. 


The  broadband  case  is  somewhat  more  difficult.  If  the  spectrum  is  perfectly  white 
(infinite  roll-off  at  the  Nyquist  frequency  and  = 0)  the  formula 


holds.  A more  general  formula  for  the  smallest  time  constant  which  holds  even  when  the 
spectrum  is  not  flat  is  given  by  (ref  4) 


T = min 

CO 


1 

2j&S(co)  fj 


where  S(co)  is  the  spectral  level  of  an  N point  FFT  of  the  reference  input. 


(A-4) 


39 


f 


More  precisely, 


S(m) 


R(j)eieM/N, 


N power 
2 cell 


(A-5) 


For  broadband  inputs,  (A-5)  may  be  rewritten. 


S(m)  = 


power 

Hz 


(A-6) 


For  a broadband  signal  which  consists  of  white  noise,  which  has  been  bandpass 
filtered  at  bandwidth  B,  the  power/Hz  of  the  spectrum  is  given  by  a^/B.  For  a sinusoid, 
the  power/cell  is  . Thus,  (A-4)  to  (A-6)  yielded. 


. N 9 

max  S(m)  =—  o~ 
m 


(A-7) 


T 


1 


+ 


2 fs  ^ 
2B 


(A-8) 


Fquation  (A-8)  is  more  general  than  ( 1 ),  and  reduces  to  ( 1 ) when  the  noise  is  white,  i.e., 
its  bandwidth  is  the  Nyquist  frequency  ^s/2. 

When  only  the  sinusoid  is  present,  (A-8)  gives  (A-2);  however,  for  broadband  noise, 
(A-3)  must  be  replaced  by 


(derived  by  substituting  (A-1 ) and  (A-8)  and  setting  = 0). 

In  the  present  experiments,  the  noise  was  low  pass  filtered  at  600  Hz  (3  dB  down- 
point).  Since  fj  = 2048,  2B  = .6  f^,  and  the  appropriate  formula  is. 


(A- 10) 
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Strictly  speaking,  several  time  constants  will  be  present  and  the  regions  of  the  output  spec- 
trum beyond  600  Hz  converge  more  slowly  (ref  4).  An  example  of  this  appears  in  figure 
A-1 . Even  after  10  time  constants,  the  higher  frequencies  (above  750  Hz)  have  not  converged. 

In  order  to  measure  the  single  time  constant  given  by  equation  (A-1 1),  the  slower 
modes  were  filtered  out  by  including  a post-filter  identical  to  the  prefilter  (fig  2).  Also, 
as  a test  of  the  above  theory,  test  1.1  was  rerun  with  the  filter  set  at  1000  Hz,  giving  a roll- 
off of  3 dB  at  860  Hz  and  compared  with  the  theoretical  formula 


The  results  are  plotted  in  figure  A-2  and,  just  as  in  figure  5,  show  excellent  agreement. 

Finally,  we  wish  to  remark  that  the  above  time  constants  were  measured  with  re- 
spect to  amplitude;  i.e.,  volts  RMS;  and  hence  are  twice  as  large  as  they  would  be  for  the 
output  power.  The  experimental  values  were  determined  by  fitting  a curve  of  the  form 
Ae“t/^  + Aq  to  the  data. 


Figure  A-1.  Output  of  ANC  with  no  post-filter  after  about  10  time  constants.  The  roll-off  of  the  prefilter 
was  3 dB  at  600  Hz. 
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APPENDIX  B - MEASUREMENT  OF  t IN  TEST  3.1 


I 


Although  we  have  been  fitting  exponentials  to  the  learning  curves,  it  is  more  con- 
venient here  to  use  the  more  exact  formula  (ref  4), 

output  = A(  1 - 1 /fjT)*' . (B-1 ) 


Then, 


20  Cn(output)  20  CnA  - n/f^r  where  f^r  > > I 

= 20  CnA  - t/r.  ^^  2) 

Then,  by  converting  the  logarithms  to  the  base  10 

20  1 

slope  of  output  (dB)  = ~ ' (B-3) 


When  more  than  two  time  constants  are  present,  the  above  technique  is  usually 
valid  provided  the  time  constants  are  widely  separated  and  the  curve  is  straight  over  a range 
of  about  6 dB.  This  may  be  seen  as  follows. 

Approximate  the  curve  by : 


A{1  -X)"  + B, 


(B-4) 


where  the  faster  time  constants  have  died  out  and  the  slower  ones  are  almost  constant 
(i.e.,  = B).  The  logarithmic  slope  of  (B-4)  is  given  by. 


nA(  1 - X)""* 

n(  1 - X)"  + B 


n 


(1  -X)  + 


B(1  -X) 
nA(l  -X)" 


(B-5) 


The  second  term  in  the  denominator  will  be  negligible,  provided  B < A(1  - X)'';  i.e., 
provided  the  magnitude  of  the  function  in  (B-4)  is  at  least  6 dB  larger  than  B.  See  figure 
30  for  an  example. 
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APPENDIX  C 


QUANTITY 


I 


I 

I 

11 


3 

1 

2 

2 

2 


EQUIPMENT  USED 
TYPE 

Universal  Adaptive  Filter,  Rockwell/Autonetics  Mod  3 

Spectrum  Analyzer,  Nicolet  Scientific  Mini-Ubiquitous 
Model  444A 

Spectrum  Analyzer,  Nicolet  Scientific  Mini-Ubiquitous 
Model  440A 

Digital  Signal  Processor,  Spectral  Dynamics  Model  SD  360 
Frequency  Translator,  Nicolet  Scientific  Model  30A/LIS 
Frequency  Translator  (Master)  Spectral  Dynamics  Model  SD  332 
Frequency  Translator  (Slave)  Spectral  Dynamics  Model  SD  332-2 
Interactive  Digital  Plotter,  Tektronics  Model  4662 
X-Y  Plotter  with  Time  Drive,  Hewlett-Packard  Model  HP  7004A 
X-Y  Plotter,  MFE  Model  715 

Pseudo-Random  Sequence  Generator,  Hewlett-Packard  Model  HPOl- 
3722A 

Random  Noise  Generator,  General  Radio,  Model  GR  1381 

Frequency  Synthesizer,  Fluke  Model  6010A 

Frequency  Synthesizer,  Rockland  Model  5100 

Frequency  Synthesizer,  Wavetek  Model  171 

Strip  Chart  Recorder,  Gultan  Model  TR  222 

True  RMS  Voltmeter,  Hewlett-Packard  Model  HP  3403C 

Pulse  Generator,  Data  Pulse  Model  101 

High  Resolution  Counter,  Hewlett-Packard  Model  HP  5307A 

Filter,  Krohn-Hite  Model  3202 

Attenuator  Set,  Hewlett-Packard  Model  350B 

Attenuator  Set,  NOSC  Manufactured 

Summing  Amplifier,  NOSC  Manufactured 

Impedance  Converter,  NOSC  Manufactured 
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